Identify Charged Higgs Boson in W ± H T associated production at LHC 
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We investigate the possibility to discover the charged Higgs via pp — > W^H^ — > I + $t + bbjj 
process at LHC, which suffers from large QCD backgrounds. We optimize the kinematic cuts to 
suppress the backgrounds, so that the reconstruction of the charged Higgs through hadronic decay 
is possible. The angular distribution of the b-jet from H decay is investigated as a way to identify 
the charged scalar from vector bosons. 
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I. INTRODUCTION 

Understanding electroweak symmetry breaking is one of the driving forces behind the undergoing experiments at the 
CERN large hadron collider (LHC). In the Standard Model (SM), the fermions and gauge bosons get masses through 
Higgs mechanism with one weak-isospin doublet Higgs field. Although SM is extremely successful in phenomenology, 
there are still remainingproblems not well understood. Extensions of SM have been considered widely. Two Higgs 
Doublet Model (2HDM)0-@ is one of the natural extensions. In this kind of models, the charged Higgs boson [H ± ) 
is of special interest, since its discovery is an unambiguous evidence for an extended Higgs sector. Therefore, the hunt 
for charged Higgs bosons plays an important role in the search for new physics at LHC. 

Currently most of the limits or constraints to the charged Higgs mass are model-dependent. The best model- 
independent direct limit from the LEP experiments is m#± > 78.6 GeV at 95% C.L.[7|, assuming only the decay 
channels H + — > cs and H + — > tv t . As the charged Higgs will contribute to flavor changing neutral currents at 
one loop level, the indirect constraints can be extracted from B-meson decays. In Type II 2HDM, the constraint is 
M H ± > 350 GeV for tan/3 larger than 1, and even stronger for smaller tan/3 Q. However, since the phases of the 
Yukawa couplings in Type III or general 2HDM are free parameters, m#± can be as low as 100 GeV0. 

At hadron colliders, the charged Higgs phenomenology has been studied widely. The main production modes are 
gb — > H~t for m H ± > m t + nib and gg — > H~tb for m H ± < m t — mj [Toj - TTsl j . The preferred decay modes are then 
H~ — > bt and H~ — > tv Ti respectively. Another interesting channel is the production in association with a W 
boson, whose leptonic decays can serve as an important trigger for the W^H^ search. This channel can also cover the 
region m H ± ~ m t . The dominant channels for W^H^ production are bb — > W^H^ at tree level and gg — > W^H^ at 
one-loop level. H /± 7J =F production at hadron colliders in Type II 2HDM and other 2HDMs has been studied in [l6rl2l| . 
It is found that, due to the large negative interference term between the triangle- and box- type quark- loop Feynman 
diagrams, the gluon-fusion cross section in MSSM is q uite small. The NLO-QCD and SUSY-QCD corrections to bb 
annihilation amplitude in MSSM are all about 10%[22L 23]. As has studied in [UHII], the hadronic decay channel of 
suffers from large QCD backgrounds, which overwhelm the charged Higgs signal over the heavy mass range that 
can be probed at LHC. In this paper, we optimize the kinematic cuts of the final state. The signal-background ratio 
is improved so that the resonance reconstruction is possible. Once the resonance of bt (tb) is detected at LHC, its 
spin should be determined. We investigate the angular distribution of the b-jet with respect to the beam direction. 
It is found that such distribution is useful to identify the charged scalar from vector bosons. 

This paper is organized as follows. In section II, the corresponding theoretical framework is briefly introduced. 
Section III is devoted to the numerical analysis of W H T production and the related SM backgrounds. In section IV, 
the angular distributions of b-jet are investigated to identify the scalar from vector bosons. Finally, a short summary 
is given. 
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(a) (b) 
FIG. 1: Feynman diagrams for W ± H^ 1 production at partonic level. 

II. THEORETICAL FRAMEWORK FOR CHARGED HIGGS 

The scalar sector of SM is not yet confirmed by experiments and it is possible to extend the Higgs structure to 
two Higgs doublets. For example, it has been shown that if one Higgs doublet is needed for the mass generation, an 
extra Higgs doublet is necessary for the Spontaneous CP violation In 2HDMs, after the spontaneous symmetry 
breaking, there remain five physical Higgs scalars, i.e., two neutral CP-even bosons ho and Ho, one neutral CP-odd 
boson A , and two charged bosons H ± . In this work we aim to study the charged Higgs phenomenology and choose 
the Type II Yukawa couplings as the working model, 
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The Yukawa couplings are related to the fermion masses, which indicates that the dominant production processes for 
the charged Higgs associated with a W boson at hadron collider are bb annihilation at tree-level and gluon fusion at 
one loop level. The bb annihilation is overwhelming for tan/3 > while large tan/3 is favored by B meson rare 

decays |25H27| . In this work we investigate the properties of charged Higgs with large tan/3 in bb annihilation process. 
The Feynman diagrams are shown in Fig. [T] The corresponding invariant amplitude square averaged over the spin 
and color of initial partons is given by 
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where a = (j>b + Pi) 2 , t = (p& — Pw) 2 1 u = (Pb ~ Ph) 2 are the Mandelstam variables and Gp is the Fermi constant. 
The Higgs propagator functions are 
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where Th .a are the higgs widths which are obtained with HDECAY[28( package. 



III. W ± H^ PRODUCTION AND CORRESPONDING BACKGROUNDS 



The total cross section for the pp — > W H T process can be written as follows 

o- = / f b (x 1 )fi(x 2 )ct(x 1 X2s)dx 1 dx2 7 
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FIG. 2: The total cross section as a function of m H ± for pp — > W ± H 11 process at (a) yfs = 7 TeV and (b) yfs = 14 TeV. 
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FIG. 3: (a) The transverse momentum distributions of the jets (ji, ]2, js, ji) with PTj x > PTj 2 > Pt J3 > Pt J4 for m H ± = 500 
GeV at yfs — 14 TeV. (b) The transverse momentum distribution of the charged lepton (solid line) and the missing transverse 
energy fir distribution (dashed line), (c) The minimal angular separation distributions between jets (solid line) and that 
between jets and the charged lepton (dashed line) . 



where ^/s is the proton- proton center of mass energy, a is the partonic level cross section of bb — > W^H^, and f q {xi) 
is the parton distribution function (PDF). In our numerical calculations, we employ CTEQ6Ll[2^ for PDF, and set 
Vtb = 1, Mw — 80.4 GeV, and nit = 172.9 GeV. In Fig. [21 the total cross sections are shown as a function of charged 
Higgs mass for tan/3 = 10, 30, and 50 at LHC. The cross section increases with tan/3. Supposing the luminosity to 
be 10 fb^ 1 at yfs = 7 TeV, one can notice that it is difficult for the charged Higgs associated with a W boson to be 
detected when its mass is above 600 GeV even for tan /3 = 50. It is easier for the charged Higgs boson to be observed 
at yfs — 14 TeV. Therefore in this work we focus on investigating the charged Higgs associated with a W boson in 
the following processes 



pp -> W~H + ->■ W~tb rvbbjj, 
pp -> W+H- -> W + tb -> l + vbbjj 



(5) 



at -y/s = 14 TeV with tan /3 = 30. 

To be more realistic, the simulation at the detector is performed by smearing the leptons and jets energies according 
to the assumption of the Gaussian resolution parametrization 
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where 5(E) / E is the energy resolution, a is a sampling term, & is a constant term, and © denotes a sum in quadrature. 
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FIG. 4: The bbjj invariant mass distributions. 



We take a = 5%, b — 0.55% for leptons and a — 100%, b = 5% for jets respectively [30j . 

The transverse momentum distributions for the four jets are shown in Fig. [3] (a). In Fig. |3] (b), the transverse 
momentum distribution for the charged lepton and the missing transverse energy ($t) distribution are displayed. In 
order to identify the isolated jet (lepton), we define the angular separation between particle i and particle j as 



A% = ^A4 + Ar??., (7) 

where Afaj = fa — faj and A 77^ = 7^ — r/j. fa (r]i) denotes the azimuthal angle (rapidity) of the related jet (or lepton). 
The corresponding distributions for AR = min(ARij) are shown in Fig.[3](c). 

The momentum of the neutrino can be reconstructed from the W mass, missing transverse momentum and neutrino 
mass, i.e., the neutrino momentum is obtained by solving the following equations 

Put = -(Pit + J2pjt), 

mw = {Pu +Pi) 2 , vl = 0, (8) 

where piT is the transverse momentum of the corresponding particle and p u (pi) is the four-momentum of neutrino 
(charged lepton). We veto the event if no solution can be found. 

In our analysis, all the hadronic jets are from the charged Higgs boson decay. As a result, the invariant mass of 
these jets can be used to reconstruct the charged Higgs boson mass. The distributions da/dMjjjj for various charged 
Higgs mass are shown in Fig. [4] 

Based on the above discussion, we employ the basic cuts (referred as cut I) 

Pit > 20 GeV, p 3 T > 30 GeV, $ T > 20 GeV, 

\rn\ < 2.5, 1^1 < 3.0, AR jm > 0.4. (9) 

For the processes Eq. (|12p with final state I + fJr + bbjj, the dominant SM backgrounds are ti, tiW, tiZ, WZjj, 
WW jj and Wjjjj, which are generated with the MadGraph 31 1 and Alpgen[32j. In the W^H^ production processes, 



four jets are from the charged Higgs decay, and three of them are from top quark decay. Therefore to purify the 
signal, one can require the invariant mass of final jets to be around the charged Higgs mass, and one top quark is 
reconstructed by three jets. Since ti is one of the predominant backgrounds, one can veto ti events if the second top 
quark can be reconstructed. Such kinds of invariant mass cut (referred as cut II) are 

I Mjjjj - m H ± I < 30 GeV, \M jj:j - m t \ < 20 GeV and \M jlv - m t \ > 20 GeV. (10) 

In order to further purify the signal, we apply a cut on the invariant mass for all of the visible particles 



M 33Jjl = J C>Zp3 + Pi) 2 > l5m H ± (11) 



together with one b-tagging (referred as cut III). 
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The cross sections for signal after each cut are listed in table |TJ We find that after all cuts, there is about 1 fb 
left for the signal process around m H ± = 500 GeV, and the backgrounds are suppressed significantly. In table HH we 
list the event numbers for the signal and background processes survived after all cuts with the integral luminosity of 
300 fb^ 1 at y/s = 14 TeV. The significance for the signal to background can reach above three sigma for m H ± > 400 
GeV. As an example, we choose m H ± = 500 GeV for the investigation in section HVl 



m H± (GeV) 


No cuts 


Cut I 


Cut I+II 


Cut I+II+III 


300 


105 


6.37 


3.33 


1.43 


400 


49.7 


4.41 


2.05 


1.11 


500 


25.7 


2.83 


1.20 


0.72 


600 


14.1 


1.81 


0.70 


0.45 


800 


4.90 


0.77 


0.26 


0.18 


1000 


1.93 


0.35 


0.10 


0.07 



TABLE I: The cross section of the signal process at \/s = 14 TeV in unit of fb. 



rn H ±(GeV) 


300 


400 


500 


600 


800 


1000 


W ± H^ 


429 


333 


216 


135 


54 


21 


tt 


14640 


3360 


696 


231 


<1 


<1 


tiW 


24 


9 


3 


2 


<1 


<1 


ttz 


30 


6 


<1 


<1 


<1 


<1 


WZjj 


69 


<1 


<1 


<1 


<1 


<1 


WWjj 


120 


117 


<1 


<1 


<1 


<1 


Wjjjj 


20940 


5970 


<1 


<1 


<1 


<1 


S/B 


0.01 


0.03 


0.31 


0.57 


>9 


>3.5 




2.27 


3.42 


8.15 


8.77 


>22 


>8.57 



TABLE II: The event numbers after all cuts for the signal and backgrounds with a integral luminosity of 300 fb 1 at yfs = 14 
TeV. The signal— background ratio and the significance are given. 



IV. IDENTIFY FROM CHARGED VECTOR BOSONS 



The reconstruction of the resonance from bt (tb) final states can straightforward lead to the conclusion that a new 
charged boson is detected. However, many theories beyond SM also predict the existence of new heavy charged vector 
bosons (e.g. W ,:lz ) which can decay to bi (tb). It can not be ignored to identify the scalar from the vector bosons with 
the identical final state. We study the following processes 

pp ->■ W~W' + ->• W~tl I'vbbjj, 

pp -> W + W'~ ->• W + tb ->■ l + vbbjj. (12) 

The corresponding transverse momentum distributions and bbjj invariant mass are similar to Fig. [3] and Fig. 01 Hence 
a proper observable has to be found to represent the spin of the new particle, which is another main aim of this paper. 
Taking into account that the scalar particle is different from the vector ones in the angular distribution of their decay 
products, we can define the angle 6 as follows 

"•-we- (13) 

where p is the 3-momentum of one of the initial proton in the laboratory frame, and is the 3-momentum of the 
b-jet which is not decay from top (anti-)quark in the H ± (or W l:t ) rest frame. The distributions ^ g related to 
H and W^- without any cuts are shown in Fig. [SI The bottom (anti-)quark in rest frame is isotropic, and has 
no correlation with the proton moving direction. However, if the new charged particle is a vector (W ,ziz ), the p£ is not 
isotropic any more. The distribution is sensitive to the chiral couplings. Such chiral couplings of W have also been 
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FIG. 5: The angular distribution of bottom quark in the (a) W process and (b) W ± W'^ process at LHC for m^ H ± w i±- ) — 
500 GeV. Solid, dashed and dotted lines respectively stand for the distributions through W'f^ , W'J 1 and W 1 ^ . The symbols 
W'^ 1 (W^ 1 ) and W't^ respectively denote the vector bosons participated in left (right) handed and pure vector type interactions. 




FIG. 6: The angular distribution of b-jet in the (a) W H^ 1 process and (b) W ± W'^ 1 process at LHC after the detector 
simulation and through all cuts, Solid, dashed and dotted lines respectively stand for the distributions through W'j~, W'^ and 
Wy* 1 . The symbols are similar to Fig. [5] 

studied in other processes [33l l34j . Fig. [6] shows the corresponding b-jet angular distributions after the smearing and 
the kinematic cuts in section iLUl One can find that the curve for iJ ± is slightly distorted by the kinematic cuts, which 
does not change the fact that the distributions corresponding to the scalar and the vector bosons are characteristically 
different. To characterize such difference, we employ the function 

f(co S e) = ^jj(i + Aco S 2 e) (i4) 

to fit the curves in Fig. [5] and [6l The values of A are listed in table Mil Obviously, A is different for scalar and vector 
bosons before or after all cuts. Therefore, the investigation of the angular distribution and the characteristic quantity 
A is helpful to discriminate the charged scalar from the vector bosons. 









W* 


w? 


No cuts 





0.69 


0.28 


0.65 


Cut I+II+III 


-0.13 


0.66 


0.14 


0.60 



TABLE III: The values of A from fitting the curves in Fig. [5] and [6] with the function Eq. {T4}. 
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V. SUMMARY 



We investigate the possibility of detecting charged Higgs production in association with W boson via pp — > 
W^HT —> I + $t + bbjj at LHC. We apply the resonance reconstruction (cut II) and resonance mass dependence 
(cut III) to suppress the QCD backgrounds. It is found that with 300 jb~ x integral luminosity at yfs = 14 TeV, 
the signal can be distinguished from the backgrounds for the charged Higgs mass around 400 GeV or larger. If a 
new resonance particle is observed, one of the key question is to identify its spin. We study the angular distributions 
for the charged Higgs and the representative vector bosons (W ,ziz ) with various chiral couplings. From the angular 
distributions and the characteristic quantity A the scalar can be distinguished from the vector bosons. The above 
analysis can be applied to discriminate the scalar from the vector bosons once the new resonance particle produced 
in association with a vector gauge boson (e.g. W or Z) is discovered at LHC. 
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